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Abstract. We present in this paper a substructure and spectroimaging study of the Coma cluster of galaxies based 
on XMM-Newton data. XMM-Newton performed a mosaic of observations of Coma to ensure a large coverage of 
the cluster. We add the different pointings together and fit elliptical beta-models to the data. We subtract the 
cluster models from the data and look for residuals, which can be interpreted as substructure. We find several 
significant structures: the well-known subgroup connected to NGC 4839 in the South- West of the cluster, and 
another substructure located between NGC 4839 and the centre of the Coma cluster. Constructing a hardness 
ratio image, which can be used as a temperature map we see that in front of this new structure the temperature 
is significantly increased (higher or equal 10 keV). We interpret this temperature enhancement as the result of 
heating as this structure falls onto the Coma cluster. We furthermore reconfirm the filament-like structure South- 
East of the cluster centre. This region is significantly cooler than the mean cluster temperature. We estimate the 
temperature of this structure to be equal or below 1 keV. A possible scenario to explain the observed features 
is stripping caused by the infall of a small group of galaxies located around the two galaxies NGC 4921 and 
NGC 4911 into the Coma cluster with a non-zero impact parameter. We also see significant X-ray depressions 
North and South-East of NGC4921, which might either be linked to tidal forces due to the merger with the 
Western structure or connected to an older cluster merger. 

Key words. Galaxies: clusters: general, intergalactic medium, general; Cosmology: miscellaneous; large-scale struc- 
ture of the Universe; X-rays: general 



"-j 1. Introduction 



The Coma cluster, also known as Abell 1656, is one of the 
best studied clusters in our Universe (see Biviano 1998 



for a historical review). It has been explored in all wave- 
lengths from radio to hard X-rays. However, despite the 
large number of observations, the physics of the Coma 
cluster and its dynamical state is not yet completely un- 
derstood. The cluster hosts a powerful radio halo (Feretti 
& Giovannini 1998) and Beppo-Sax data recently revealed 



the existence of hard non-thermal X-ray emission (Fusco- 
Femiano et al. 



1999) 



Extensive work has been carried out to measure the ve- 
locity dispersion and distribution of the galaxies located 
in the Coma cluster. CoUess & Dunn (1996) for example 



published a catalog of 552 redshifts and determined the 
mean redshift of the cluster to be cz = 6853 km/sec and a 
veloci ty dispersion of cTcz = 1082 km/sec. CoUess & Dunn 
( |l996|) find and confirm the presence of substructure in the 



galaxy distribution: there is a sub-group located around 
NGC 4839, which lies to the South- West of the cluster. 
The existence of this s ubstructure was originally foun d 
by Mellier et al. ( |l988|) and Merritt & Trimbley ( |l994|) 
Furthermore, it was noted by Fitchett & Webster (1987) 
that the two central galaxies, NGC 4889 and NGC 4874 
have a significant difference in velocity, which is indicative 
of a recent merger. It is probable that NGC 4889 belonged 
to a subgroup, which recently merged with the main clus- 
ter and that NGC 4874 is the central dominant galaxy of 
the main cluster. 

There are a number of papers addressing the thermal 
X-ray emission of the Coma cluster which comes from the 
hot thermal intracluster medium (ICM). White, Briel & 
Henry (199S) discussed substructures in the Coma cluster 
based on ROSAT-data. Briel, Henry & Bohringer (|1992D 
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analyzed the ROSAT All-Sky-Survey data and determined 
the X-ray surface brightness profile of the cluster out to 
a radius of roughly 100 arcmin. Hughes ( p^89| ) performed 
an extensive mass analysis on the Coma cluster based on 
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X-ray and optical data. Later work by Vikhlinin, Forman 
& Jones ( |1997 ) revealed the existence of a filament- like 
substructure, which is most likely linked to NGC 4911, 
a bright spiral galaxy located South-East with respect to 
the centre of the cluster. Spectroimaging studies based 
on ASCA-data constrained the X-ray temperature distri- 
bution of th e ICM. Honda et al. ( |1997| ) and Watanabe 
et al. (199£) presented temperature maps of the Coma 
cluster based on ASCA-data, which clearly showed that 
the ICM is not isothermal. Furthermore, Donnelly et al. 
( p^99| ) saw indications for a hot spot North of the centre of 
the Coma cluster. Recently, first results of the Coma clus- 
ter based on XMM-Newton data were published: Briel et 
al. ( 2001 ) discussed the overall morphology of the clus- 
ter and showed the existence of several point-like sources, 
which are linked to individual galaxies in the cluster and 
which have not been accounted for in previous studies. 
Arnaud et al. ( 2001a ) presented a temperature map of the 
central part of the cluster based on XMM-Newton data 
with unprecedented spatial resolution. This map showed 
again variations of the X-ray temperature (see Jansen et 
al. pool] for an overview of XMM-Newton and Turner et 
and Strueder et al. 



al. 2001 



2001 for an overview of the 



EPIC detectors). However, the previously found hot spot 
found by Donnelly et al. 1999 was not confirmed. Lastly, 
Neumann et al. ( 200l| ) presented strong evidence based on 
XMM-Newton data that the substructure located around 
NGC 4839 (see above) is on its first infall onto the cluster 
from the South- West and has not yet passed its centre. 
This is in contradiction with previous work by Burns et 
al. 1994, who suggested that the group around NGC 4839 



has already passed the centre of the Coma cluster once. 

The aim of this paper is to perform a substructure and 
spectroimaging study of the Coma cluster based on XMM- 
Newton data, to better understand the underlying physics. 
The work presented here is our first attempt to discuss 
the overall dynamics of the cluster, and is by no means 
exhaustive. It aims to show the large potential provided 
by the data obtained with XMM-Newton . The paper is 
organized as follows: after the introduction we present the 
observations in Sec. 2. This is followed by the description of 
our data treatment in Sec. 3. In Sec. 4 we present our spa- 
tial analysis and in Sec. 5 we describe our spectro-imaging 
results. In Sec. 6 we discuss our results and finish in Sec. 7 
with the Conclusion. 

Throughout this paper, we assume Hq = 50 km s^^ 
Mpc-\ = 1 and go = 1/2. At a redshift of z = 0.0232, 
1' corresponds to 38.9 kpc. The cluster is sufficiently close 
so that variations of the geometry of the universe do not 
matter. 



2. Observations 

The size of the Coma cluster (roughly 3 degrees in di- 
ameter) greatly exceeds the field-of-view (FOV) of XMM- 
Newton . Therefore the observations of Coma were per- 
formed in mosaic mode with several exposures pointing at 
different regions of the cluster. The summary of all obser- 



vations is shown in Tab.^. The exposure times indicated 
are the effective observing times after flare rejection. We 
apply the same method and thresholds for flare rejection 
as described in Majerowicz, Neumann & Reiprich (2002) 
for the pn and MOS-cameras. During the observations the 
pn-camera was mostly operating in extended-full-frame 
mode. There was one additional Coma pointing, in which 
only the pn-camera was active (Coma 2), for which we 
were not able to run the pipeline processing at the time 
of the analysis. Thus we discard this observation in the 
following. 

3. Data treatment 

Our analysis presented here is based on SAS-version 5.1 
and concerns only events in the field-of-view (FOV). 



3.1. Vignetting correction 

In order to correct for vignetting we apply the photon 
weighting method described in Arnaud et al. (2001b). 
To outline it briefly: to each observed event at energy 
E falling onto the detector position x,y we attribute a 
correcting factor: 



UJ 



effective area(a; = 0,y = 0, E) 
effective area(a;, y, E) 



(1) 



the effective area takes into account the mirror surface, 
the RGA structure (in the case of the MOS-cameras), the 
filter transmission as well as the quantum efficiency of 
the detectors. We assume in our anaylsis that the filter 
transmission and the quantum efficiency are flat over the 
entire field of the detector, x ~ Q,y — stand for on-axis 
position. 

All images presented in this study are vignetting cor- 
rected. For this we create images by adding the different 
weight factors of the photons into image pixels. The sta- 
tistical uncertainty per image pixel is calculated as 



a = ^Sr^iu;2 (2) 
where n is the number of photons per pixel. 

3.2. Merging of the different exposures 

While merging the data from different exposures we ob- 
served certain offsets of the order of a few arcseconds for 
which we corrected by looking at point sources, which were 
visible in two or more pointings. We shifted the pointings 
so that the same point source, visible in different expo- 
sures, appeared in each pointing at the same sky coor- 
dinate position. We estimate the remaining errors in the 
order of a few arcseconds. 

Fig.|l| shows the resulting countrate image in the en- 
ergy band 0.5-2.0 kcV after merging the different observa- 
tions together (MOS and pn) and subtracting the back- 
ground. For the background we used the Blank-Sky ob- 
servations provided by D. Lumb and retrievable from the 
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pointing 


camera 


exposure 


RA 


Dec. 


date 


no. 




in sec 


in deg. 


in deg. 


mm/yy 


10 


ml 


18576 


194.910 


28.1278 


06/00 




m2 


18584 










pn 


17500 








11 


ml 


13845 


194.652 


28.3989 


06/00 




m2 


15237 










pn 


13754 








1 


ml 


26559 


194.199 


27.4019 


06/00 




m2 


26594 










pn 


23732 








2 


ml 


26046 


194.427 


27.7272 


12/00 




m2 


26062 










pn 


15000 








3 


ml 


12519 


194.634 


27.4033 


06/00 




m2 


11523 










pn 


11730 








4 


ml 


476 


195.019 


27.5233 


06/00 




m2 


582 










pn 


2970 








5 


ml 


14977 


194.865 


27.7814 


05/00 




m2 


16485 










pn 


14211 








611 


ml 


14943 


194.708 


27.9811 


12/00 




m2 


14942 










pn 


8874 








6 


ml 


6548 


194.708 


27.9811 


06/00 




m2 


6744 










pn 


3969 








711 


ml 


23755 


194.365 


28.1447 


12/00 




m2 


23245 










pn 


8874 








7 


ml 


3153 


194.365 


28.1447 


06/00 




m2 


2960 










pn 


358 








8 


ml 


23455 


195.357 


27.7314 


12/00 




m2 


24172 










pn 


18430 








9 


ml 


19951 


195.136 


27.9497 


06/00 




m2 


20361 










pn 


15443 








bgd 


ml 


21145 


195.404 


27.3311 


06/00 




m2 


21563 










pn 


17534 








can 


ml 


16103 


194.945 


27.9500 


05/00 




m2 


15894 










pn 


14595 









Table 1. The summary of exposure times and coordinates 
for the different Coma pointings. 



XMM-Newton web-pages. The background subtraction is 
not important for the construction of the countrate im- 
ages, however, it is important for the determination of the 
hardness ratio map (see below). 

In order to obtain Fig.|^ we perform the following data 
treatment: 

— Merge the photon weight images of the individual MOS 
(M0S1-I-M0S2) pointings together (image: mergemos) 
in the energy range 0.5-2.0 keV 



Create the corresponding MOS mosaic exposure map 
(image: mergemosexpo) , which consists of the sum of 
the different detector masks|] per observation and cam- 
era (MOSl and M0S2), which are multiplied by the 
observing time of the corresponding exposure before 
summing. 

Divide mergemos by mergemosexpo = countratemos 
Create corresponding background countrate image for 
MOS (image: hackgroundmos) in the same energy 
range as in mergemos. For each MOS exposure we con- 
struct a vignetting corrected background exposure, for 
which we use as input the blank sky observations com- 
piled by D. Lumb. The corresponding photon weight 
images are added together. We construct the corre- 
sponding exposure map of the background and divide 
the photon weight image by the exposure map. The re- 
sult is a background countrate image hackgroundmos. 
countratemos - hackgroundmos = sourcemos 
Merge the photon weight images of the pn images (im- 
age: mergepn) in the same energy range as was used 
for MOS. 

Create the corresponding pn exposure map (image: 
mergepnexpo) in the same way as done for MOS 
Create background countrate image for pn in the same 
manner as for MOS (image: backgroundpn) 
Divide mergepn by mergepnexpo — countratepn 
countratepn - backgroundpn = sourcepn 
sourcemos + sourcepn = total background- subtracted 
countrate image = Figj^ 



For all observations MOS and pn were operated at the 
same time and had, perhaps with the exception of pointing 
no. 4, sufficiently long exposure times to give for each type 
of instrument an adequate countrate estimate per pixel. 
Therefore it is correct to add together the different coun- 
trate images obtained by the different detectors. If for one 
observation only one type of camera had been available 
we could not add the countrate images since this individ- 
ual exposure would create biases in our mosaic image due 
to the different sensitivities and spectral responses of the 
detector types. 

We want to stress that the Coma pointings were per- 
formed in medium filter mode while the data used for the 
background are compiled from observations mainly using 
thin filter. However, since we only use photons at energies 
above 0.5 keV we expect the difference of the different fil- 
ters to be negligible. Furthermore, at regions with high 
signal-to-noise ratio we expect that the background sub- 
traction is not important. 



^ The detector mask represents the geometry of the detector 
which includes gaps between the CCD's, the region of FOV, 
and excludes bad pixels, which were identified in the pipeline 
processing 
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Fig. 1. The Coma cluster in the 0.5-2.0 keV band. The image is the sum of aU the different exposures (see Tab. I) of 
the different MOS and pn cameras (see text). The image is background and vignetting corrected. We apphed a small 
median filter to mask out effects caused by detector gaps. Colour scale: dark blue regions correspond to countrates of 
0.016 cts/sec/arcmin^ and white regions to countrates > 0.288cts/sec/arcmin^. 



4. Spatial analysis 

4.1. Substructure analysis 

Our aim in this study is to look for substructure and tem- 
perature variations in order to better understand the dy- 
namical state of the Coma cluster. 

For measuring substructure we fit an elliptical beta- 
model to the cluster images, which describes a relaxed 
cluster without substructure . Subsequently, we create 
a synthetic cluster image based on the best fit parame- 
ters of the beta-model. This model is subtracted from the 



data. The residuals remaining after subtraction, which in- 
dicate substructure, are examined for their significance. 
This method has been applied to useful effect in sev eral 
other studies, e.g., Neumann & Bohringer 1997 199£ and 



Neumann 199E 



In the following we describe the different steps of our 
substructure analysis in detail. 
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MOS 


pn 


observed background 


1.2 ±0.3 


2.9 ± 1.1 


countrate in 10~*counts/sec/pixel 






fitted background 


1.2 


3.8 


countrate in 10~* counts/sec/pixel 







Table 2. The observed and fitted background values of 
the Coma exposures. The observed background is deter- 
mined from the blank-sky- field observations. The fitted 
background value is the result from the elliptical beta- 
model fitting. The error bars are 68% confidence level. 



4.1.1. Beta-model 

In fitting elliptical beta-models to the XMM-Newton data, 
we take as input the vignetting corrected images in the 0.5- 
2.0 keV energy band. This energy range was chosen since it 
provides the best signal-to-noise ratios (see Scharf 20021 ). 
For this analysis we did not subtract the background, since 
we include it as fit parameter. We compare (see below and 
Tab.||) the fitted background values with the countrate 
values observed in the background observations. This is a 
consistency check for the applied background subtraction. 

We encounter small number Poisson statistics in the 
pixels located at the outskirts of the cluster. This pro- 
duces non-symmetric error bars for the pixel values. Since 
we use statistics for our beta-model fitting, which as- 
sumes a symmetric Gaussian probability distribution for 
the photons, we have to correct for the asymmetry of the 
error bars. An easy way to do this is to apply a Gaussian 
filter (cr = 15") to the image before fitting the model. This 
procedure has proven to be valid and gives good results 
as shown in previous studies (e.g. Neumann & Bohringer 



1997 , 1999 ; Neumann 1999 ) . In order to account properly 
for the vignetting, in our calculation of the error bars per 
pixel we quadratically add the different weighting factors 
for each event in each pixel before Gaussian filtering (see 
eq(2)), and use this image as estimate for the uncertainties 
in each pixel. 

To check for possible instrument-dependent variations 
we treat the MOS- and the pn-data separately. The fit 
results are displayed in Tab.^. They are in good agree- 
ment with each other and with previous studies, like Briel, 
Henry & Bohringer (1992) . Since we are only interested in 



the best-fit model in order to subtract it from the cluster 
data and check for residuals, we do not determine error 
bars on the fit parameters, which would be time consum- 
ing since the error bars would have to be estimated via a 
Monte-Carlo approach. 

The centre of both the MOS and the pn-data fit coin- 
cide relatively well with the position of the central galaxy 
NGC 4874, at: RA=194.898 and Dec.=27.963. The dis- 
tance between the model centre and this galaxy is less 
than 2 arcmin. 



4.1.2. Significance of residual map after beta-model 
subtraction 



We subtract the cluster model based on the best fit pa- 
rameters of the beta-model and calculate the correspond- 
ing significance of the residuals for each instrument type 
independently. The exact method we use for calculating 
the significances is described in Neumann & Bohringer 
(1997). For our residual calculation we apply a Gaussian 
filter with a — 1' . We obtain two significance maps, one for 
the MOS-cameras, and the other for the pn-camera. The 
residuals obtained for each type of camera are very similar. 
We quadratically add the two different significance maps 
in order to obtain a residual image with better statistics. 
The result is shown in Figj^. The image over which we 
display the significance contours is the same as in Fig.|l|. 

We want to stress that the resulting residuals with 
a significance equal or higher than -1-10 cr are extremely 
stable against changing the parameters of the subtraced 
model. 



4.1.3. Results 

The residual map shown in Fig.H shows several significant 
structures: 

— the extended X-ray emission around and behind 
NGC 4839, which was already detected and described 
in previous papers. 

— a significant (the lowest contour is 5a) structure in 
the West, which is located between the Coma centre 
and the subgroup around NGC 4839 at a distance of 
about 0.9 Mpc. The structure is elongated in North- 
South direction with an extent of about 50 arcmin, 
which corresponds to about 1.7 Mpc. The width of the 
structure in East- West direction is about 20 arcmin or 
0.6 Mpc. The structure shows two extended maxima 
at RA=12:58:1.26 and Dec.=27:44:32.25 for the first 
maximum and RA=12:58:16.88 and Dec.=28:6:17.68 
for the second. The entire structure is slightly bent 
tangentially to the Coma centre. 

— a residual with a filamentary shape is found to the 
East of the cluster centre. A similar residual at the 
sam e loc ation was already presented by Vikhlinin et 
al. (|1997| ) in a wavelet analysis of ROSAT data. The 
structure is elongated in the East- West direction and 
bends to the North at its Eastern end. The projected 
length of this filamenary structure is in the order of 
28 arcminutes (which corresponds to roughly 1 Mpc) 
and has a width in the order of 6 arcminutes (about 
200 kpc). 

— a deficit in X-ray emission is detected North and 
South-East of the filamentary Eastern structure (see 
above). The precise significance and shape of these 
X-ray decrements is hard to determine (it is around 
-20cr), since it is very sensitive to the adopted back- 
ground level, however, the signal is persistent. This 
deficit in X-ray emission can already be observed in 
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camera 


So 


Xa 


yo 






/3 




cts s^^ arcmin~^ 


in deg. (J2000) 


in deg. (J2000) 


in arcmin 


in arcmin 




MOS 


0.094 


194.929 


27.9586 


11.1 


9.0 


0.72 


pn 


0.28 


194.927 


27.9597 


12.1 


11.2 


0.85 



Table 3. The parameters from the elliptical beta-model fit. 5*0 is the central intensity, xq and yo are the central 
coordinates in RA. and Dec. respectively. Vcmaj and rcmin are the two core radii in major and minor axis direction. 



Fig.|l| comparing the overall emission South-East of the 
cluster with the overall emission North and West of 
the Coma cluster. We do not show contours of these 
negative residuals in Fig.^ since their shape strongly 
depends on the adopted background level. 
— an excess of emission in the centre, which is linked 
to the two dominant central galaxies NGC 4874 and 
NGC 4889. These residuals are signs for the two galax- 
ies to have deep individual gravitational potentials, 
which create small gas halos within the larger scale 
of the cluster. 

One important question is whether the detected resid- 
uals might not be an artefact due to improper modelling, 
for example the residuals could be due to the wrong de- 
termination of position angle of the subtracted model. In 
order to investigate this possibility we construct cluster 
models with different position angles which are subse- 
quently subtracted from the Coma data. Even when we 
change the position angle by 90 degrees, the structures 
with positive residuals keep their overall morphology and 
significance. Since the centre of the model is very close to 
the central galaxy NGC 4874, it is most likely that the 
residuals cannot be caused by incorrect positioning of the 
centre of the subtracted model. 



4.2. Spectro-imaging analysis 

In order to measure the temperature distribution of the 
ICM we determine the hardness ratio map for the Coma 
cluster. We make use of Fig.0 and of the countrate image 
in the energy band 2-5 keV. 

Both o f th ese images are constructed as described 
in Sect. 3.2 and are fully background-subtracted. 
Background subtraction is crucial to the construction 
of the hardness ratio map; without this correction 



background-dependent variations would appear in the fi- 
nal image. 

In order to avoid artefacts due to photon statistics, 
we apply a median filter to the images in the 0.5-2.0 keV 
and the 2.0-5.0 keV band. We filter over a region of 12 
by 12 pixels (pixel size is 15" x 15") before we divide the 
image obtained in the higher energy band by the image 
in the low energy band. Applying a median filter not only 
helps for the statistics but also filters out the emission of 
individual sources such as galaxies or AGN. Since there 
are insufficient statistics in the outskirts of the cluster to 
reliably constrain the hardness ratio, we restrict ourselves 
to regions in which the countrate in the energy band 2- 
5 keV band is above 2 x 10~^ cts/sec/pixel. This is a very 



conservative limit (see below) which helps avoid artefacts 
linked to imperfect background subtraction or to small 
number statistics. 



4.2.1. Results 

The hardness ratio image is shown in in Figure |. The 
image is in three colours; red shows temperatures below 
8 keV, yellow regions between 8-10 keV and white are 
regions with temperatures above 10 keV. Our results are 
in good agreement with those obtained by Arnaud et al. 



( ^OOlaD and Briel et al. ( ^00 1] ). 

Are the differences in hardness ratio significant or not? 
In order to evaluate the errors, we perform the following 
approximation: 

we only look at regions in the hardness ratio image in 
which the countrate in the 2-5 keV band is larger or equal 
to 2 X 10"'' cts/sec per pixel (the pixel size is 15" x 15"). 
Assuming a mean exposure of 20 ksec per camera, we ob- 
tain 8 photons per pixel]^ . Since the observed mean hard- 
ness ratio is roughly 0.3 (which corresponds to a plasma 
temperature of about 8 keV), the countrate in the 0.5- 
2 keV band is a factor of about 3 higher than the coun- 
trate in the higher energy band. Thus, we obtain roughly 
24 photons per pixel in a 20 ksec exposure. This calcu- 
lation takes only into account the source countrate. For 
the proper error calculation we also have to include the 
background countrate which we subtracted via the blank 
sky observations. The total background countrate in the 
0.5-2.0 keV band is 4.1 x 10~*cts/sec/pixel and for the 
2.0-5.0 keV band 2.3 x 10~^cts/sec/pixel. Furthermore, 
for the hardness ratio image we applied a median filter 
derived from the pixel values in a window of 25 x 25 im- 
age pixels. If we assume that the median filter has the 
same error as if one averages over all pixels, a good ap- 
proximation, we can determine the errors by calculating 
the total number of photons and using error propagation 
including the background. We obtain a 1 — cr-error of 0.01 
in the hardness ratio. This corresponds to about 0.5 keV 
for temperatures below 10 keV. We want to stress that 
this is the error estimate for the outer region. In the inner 
parts, where the statistics are much better, the errors are 
much smaller. As an example: from the outer parts from 
the cluster, to the centre, the increase of X-ray intensity 
is of the order of a factor of 20. 



^ this takes into account the fact that we calculated the coun- 
trate assuming only one MOS camera with a total exposure 
time of MOSl plus MOS2. 
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Fig. 2. The image is the same as in Fig.|i| only with different image value cuts and a linear look-up-table (lut) instead of a 
logarithmic lut. The contours are the significances of the residuals above main relaxed cluster and background emission 
(see also text). The lowest contour and the step width between two contours are each 5 a. Colour scale: darkblue regions 
correspond to a countrate of 0.0192 cts/sec/arcmin^ and white regions to countrates > 0.32 cts/sec/arcmin^. 



One final point concerning the results: the hardness 
ratio image shows two "hot spots" to the South. The sig- 
nificance of these two points is very low since they lie both 
in pointing no. 4, which has an exceptionally low exposure 
time (see Tab.|l|). It is thus quite likely that these features 
are statistical artefacts. 



5. Discussion 

5.1. Infall of a substructure located between Coma 
center and NGC 4839 

In Figure ||, there is an obvious, large residual to the West 
of the cluster centre. The region lying between this struc- 
ture and the Coma centre shows a temperature increase to 
temperatures equal to or above 10 keV. This suggests that 
it is heated via compression or via shock waves, which were 
created during the infall of this sub-structure onto the 
Coma centre. The creation of shock waves is observed in 
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Fig. 3. The hardness ratio image obtained by dividing the image in the 2-5 keV band by the image in the 0.5-2 keV band 
(see also text). Assuming an isothermal spectrum, red regions correspond to parts of the cluster with temperatures 
below 8 keV, yellow regions correspond to kT >8 keV and white regions have temperatures above 10 keV. The contours 
are the significance of the residuals, see also Fig.^. 



hydrodynamic simulations of gal axy c luster mergers (see 
for example Schindler & Miiller 1993 ). It seems unlikely 
that the structure is a chance alignment in the line-of- 
sight. 

The two maxima in this western structure could be 
caused by a disruption of a galaxy group during its infall 
onto the cluster centre into several smaller units. It is rela- 
tively unlikely that two galaxy groups are falling onto the 
Coma cluster at exactly the same time. 

Previous work has shown that the galaxy group around 
NGC 4839 is falling onto the Coma cluster (see for ex- 
ample Neumann et al. 2001 and references therein). The 
infall direction coincides with a filament which connects 
the Coma cluster with Abell 1367, and which is part of 
the 'Great Wall'. One may ask if the structure between 



NGC 4839 and the cluster centre is also falling in from 
the same direction. If a subcluster infalls from this direc- 
tion, it will heat up the gas between it and the main cluster 
via compression. This will lead to an increase in the ram 
pressure encountered by the subcluster. Pressure increases 
linearly with temperature for an ideal gas (P = nkT) : the 
relative difference in temperature between the hot region 
and the main cluster is ~ 20%, indicating a similar in- 
crease in pressure, and which may be sufficient to slow 
down the infalling gas. It is possible that, due to inertia, 
the gas in the subcluster is not abruptly slowed down, 
but instead may move in a direction where there is no 
increase in pressure. This movement could be tangential 
from South to North. In this case, the hot region to the 



D.M. Neumann et al.: The dynamical state of the Coma cluster 



9 



North can naturally be explained as due to compression 
of the remaining gas in the subcluster. 

The increase of temperature towards the North has 
been seen i n ASC A-data (see Watanabe et al. 1999 and 
Honda et al. 1996). It was interpreted as being due to heat- 
ing by a shock wave created from a previous merger from 
a galaxy group located around NGC 4889, which is now 
one of the two main galaxies in the centre of the Coma 
cluster. It is not possible at present to distinguish between 
the two explanations. 

Examination of DSS plates does not show any obvi- 
ous overdensity of galaxies coinciding with the Western 
residual structure. In this case, one might ask whether 
this structure is the remaining gas from the subcluster 
which was once associated with NGC 4889. Ram pres- 
sure is known to act much more efficiently on gas than on 
galaxies, which would provide a natural explanation for 
the differences in location of the gas and the galaxy. The 
effectiveness of ram pressure stripping can be seen in the 
case of the subgroup around NGC 4839, in which the gas is 
clearly lagging b ehind the galaxies on the infall trajectory 
(Neumann et al. 2001 ). Taking the projected distance be- 
tween NGC 4889 and the western substructure (~ 1 Mpc) 
and a typical cluster merger or sound crossing timescale of 
10^ years, we can calculate the mean difference in velocity 
between the galaxy and the gas. This value is close to the 
velocity dispersion of Coma (in the order of 1000 km/sec 
- Colless & Dunn 1996), which is not surprising given that 
typical cluster sizes are of order of a few Mpc. The fea- 
sibility of this scenario can only be assessed with high 
resolution hydrodynamic simulations, which take into ac- 
count the differences of ram pressure between the galaxies 
and the gas of the merging subgroup. 

5.2. The residual to the East 

This compact region clearly shows lower values in the 
hardness ratio image, indicating that this part of the clus- 
ter is cooler than the mean cluster temperature. 

The filamentary form of the structure suggests that it 
may be linked to some sort of infall process in connection 
with ram pressure stripping. Since ram pressure stripping 
works efficiently only in the centre of the cluster, where 
the gas density is high, we suppose that this structure is 
relatively close to the centre of the cluster and not a pro- 
jection effect. NGC 4911, a bright galaxy which lies in the 
region of the substructure, shows a relative difference in 
velocity with respect to the Coma cluster of 1000 km/sec, 
which is close to the velocity dispersion of the cluster as 
a whole. 



to estimate the mass of gas of this structure. A high value 
for the mass would suggest that the gas comes from an 
intragroup medium, a low value would indicate that the 
origin of the gas is from gas located inside the galaxies. 

The structure has a projected distance from the Coma 
centre of roughly 1 Mpc and enhances the observed fiux 
with respect to the fitted beta-model by about 30%. Using 
a length of 1 Mpc, a diameter of 200 kpc, and assuming 
cylindrical symmetry (thus a thickness of the structure in 
the line-of-sight of 200 kpc), the volume of the substruc- 
ture Vs - 0.04 Mpc^ 

In order to calculate the gas mass of this structure 
we have to determine the contributions of the main clus- 
ter gas and the gas in the substructure. The contribu- 
tion of the main cluster gas has to be integrated along 
the line-of-sight. For this we make use of the previously 
determined betamodel parameters by Briel et al. ( 1992 ) 
(/3 = 0.75, rc = 0.42 Mpc, = 2.9 x IQ-^cm-^) . The 
determination of the gas mass of the substructure depends 
on its location with respect to the cluster. If located at a 
large distance to the cluster we can simply add the two 
X-ray contibutions of main cluster and substructure. In 
this case we find an electron density of the subgroup of 
Ties — 1.1 X 10^^ cm^'^. Integrating this density over the 
previously calculated volume we find a total gas mass for 



the structure of M„ 



8 X 10 A/q. If we assume that 



the structure is very close to the centre and that in fact the 
physical distance towards the Coma centre is close to its 
projected distance we need to take into account that the 
gas in the substructure is at the same position than the 
gas of the main cluster, which enhances the X-ray emis- 
sion in this region^ In this case the required gas density 
for the subgroup is lower than in the case of no physi- 
cal connection between the two components. We obtain 
in this case an electron density of ries = 7.5 x 10~'*cm~^, 



5 X 10 Combining the 



and a gas mass of Mgas. 
two cases as the two possible extremes we find that the 
gas mass of the substructure lies between 5 x 10^^ -^0< 
Mgas,s < 8 X 10^^ -^^0- This value for the gas mass is at 
the low end for a group of galaxies (Mulchaey 2000) and 
suggests that a part of the gas might also come from the 
interstellar medium of galaxies themselves. The obtained 
values are in very good agreement with the value obtained 
by Vikhlinen et al. |1997| with 5 x 1O"M0. 

If we assume that the gas mass of the substructure 
represents roughly ten percent of the total mass of the 
substructure, which is of the order of magnitude of what 
is typ ically found (see for example David, Jones & Forman 
1995 ), then we obtain for the total mass Mtot,s ~ 5 — 8 x 
10^^ -^^0- This is a few times larger than the mass of 



5.2.1. Gas mass estimate of substructure 

To determine the origin of the filamentary Eastern struc- 
ture, and more specifically to know whether the the ob- 
served stripped gas comes from an intragroup medium or 
from an interstellar medium linked to galaxies, we attempt 



^ If the two regions are only a chance alignment in the line- 
of-sight the X-ray emission (xem) in this region follows: xem oc 
n^s + n^m where riem is the electron density of the main cluster 
gas. In the case that the two regions are located at the same 
place xem behaves as: xem oc rietot = {ries + nem)'^ ~ n^s + 
2neBnem + «-em > '^L + . The factor 2nesnem boosts the 
X-ray emission if the two regions are physically connected. 
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a massive galaxy and again somewhat at the low end of 
the mass of groups of ga laxies with typically a few 10^^ 
M0(e.g. Mulchaey pOO0| ). NGC 4911 and NGC 4921 lie 
in this sub-structure and are by far the brightest galaxies 
(see Fig.^ in this region. Dow & White 1995 give a value 
for the optical luminosity of NGC4911 of logLsii'a Lq) 
= 10.89. Since NGC 4921 has a similar magnitude (see Doi 
et al. 1995), we can assume a similar luminosity. Assuming 
a mass-to-optical light ratio ^ 10, we obtain a total mass 
of NGC 4911 and NGC 4921 of ~ lO^^ Mq, which is of 
the same order of magnitude as the value required for the 
entire filamentary substructure. This indicates that the 
substructure is largely dominated by the two galaxies. 



Matsushita ( 2001 ) recently presented a correlation be- 
tween mass of the hot interstellar medium and optical 
luminosity Lb for late type galaxies. Applying this cor- 
relation to NGC 4921 and NGC 4911 we find a mass for 
the hot phase in the order of 10^° ^0- This is smaller 
than the value of a few times 10^^ -^0- If correct, this 
implies that a large fraction of the gas comes indeed from 
an intra-group medium. 

The X-ray structure we see likely originates from gas 
in the surrounding intra-group medium of a small group 
of galaxies plus, to some fraction, which we expect to 
be in ther order of a few up to 10 or 20 per cent, 
from the stripped interstellar medium of the two galax- 
ies, NGC 4911 and NGC 4921. 

A cursory examination of Figure || would suggest that 
the principal direction of motion is to the West with per- 
haps an additional component to the North. Figure ^ 
shows the X-ray residual over the optical image from the 
DSS, indicating that NGC 4921 is at the Eastern end of 
the structure and NGC 4911 is right at the position where 
the structure bends. The only way we can explain this 
alignment is by assuming a direction of motion to the East, 
with NGC 4921 at the head of the structure and NGC 4911 
behind. The different directions of the stripped tails may 
indicate different velocity directions for the galaxies, with 
NGC 4911 only having a projected velocity component 
towards the East while NGC 4921 has an additional pro- 
jected component to the North. This could indicate that 
the structure in which NGC 4911 and NGC 4921 are 
bound together changes its direction of velocity with time. 

Are these two galaxies bound together? If they are, 
they should have similar measured velocities or redshifts. 
A combination of NASA Extragalactic Database (NED) 
and published data (Smith et al. 200C ) indicates velocities 
of cz = 7973 kni-i and cz = 7560 km^^ for NGC 4911 and 
NGC 4921, respectively, suggesting that they are indeed 
gravitationally bound. (Note that the NED velocity for 
NGC 4921 is incorrect, having been superseded by the 
Smith et al. measurement.) Since the two galaxies have 
larger measured velocities than the Coma cluster itself, 
it is reasonable to assume that the structure is located 
between us and the Coma cluster, thus lying somewhat 
in front of the Coma centre (although there is the small 
possibility that they are background). 



We can estimate the temperature of the sub-structure 
by assuming that it contributes about 30% to the total 
flux in the region, where it is located (see above). We per- 
form simulations of spectra, assuming a main cluster tem- 
perature of 8 keV and a sub-group temperature of 1 keV 
and calculate the corresponding hardness ratio. We find 
a hardness ratio of 0.29, in agreement with the hardness 
ratio values in the region of the substructure. If interac- 
tions occur and are important for the emission mechanism, 
we expect the cold phase to be at a temperature close to 
or below 1 keV. Thus, 1 keV seems to be a good upper 
limit for the temperature of the structure. Temperatures 
around or below 1 keV were found to be typical values for 
compact grou ps of galaxies (see for example Helsdon & 
Ponman pOOOl) . 

In the scenario described above, the two galaxies must 
have a tangential component with respect to the cen- 
tre of the main cluster, which implies that the merger 
is not a simple two-body encounter. The Coma cluster 
is surrounded by fi lamentary structures (see for example 
Huchra et al. 1990), and their gravitational pull could ex- 
plain the non-zero impact parameter of the two galaxies 
on their infall onto the cluster. 



5.3. The two dominant central galaxies 



Vikhlinin et al. (2001) have shown based on a study of 
Chandra data that the temperatures of the two dominant 
central galaxies NGC 4874 and NGC 4889 are significantly 
cooler than the surrounding gas. In our study we do not 
see a temperature drop in these regions, primarily because 
we average over a too wide area in our analysis to resolve 
this cooler gas. 

5.4. The X-ray deficit in the South and East of 
NGC 4911 and NGC 4921 



The X-ray deficit in the South and East is difficult to ex- 
plain. Comparing the overall morphology of the gas and 
temperature distribution with recent studies based on hy- 
drodynamical simulatio ns, su ch as Roettiger et al. ( 1997 ) 
and Ricker & Sarazin (2001), we cannot find a scenario 
which matches all observed features. We therefore think 
that this decrement in X-ray emission is either due to tidal 
forces linked to the merger with the Western structure 0, 
or to an older merger, possible connected to the Northern 
part of the Western substructure. This might be the rem- 
nant of a sub group, which already passed the core once. 
However, the physics of this deficit in emission might be 
more complex. We want to stress, again, that this decre- 
ment is not an artefact of our beta-model fitting. Changing 
the parameters for the model (position angle, core radius. 



^ in this case the gas could follow the mass enhancement 
linked to the structure falling onto the centre from the West 
and provoking a diminuition of gas density in the East and 
South East 



Fig. 4. The contours of the residuals superposed on the image of the DSS. The contours are the same as in Fig.g and 
Fig.|. 



centre position and background) over a considerable range 
does not alter the result (see above). 



6. Conclusion 

We performed a substructure and spectroimaging study 
based on the XMM-Newton data of the Coma cluster. The 
data show the existence of substructures which are very 
likely interacting with the main cluster. The combination 
of spatial and spectral information based on data obtained 
with XMM-Newton allows us to put global constraints on 
the dynamical state of the Coma cluster. These are impor- 
tant for our understanding of the physics of clusters. The 
results presented here show the large potential of future 
studies of clusters of galaxies observed with XMM-Newton 
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